Excitations from a Bose-Einstein condensate of magnons in coupled spin ladders. 
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The weakly coupled quasi-one-dimensional spin ladder compound (CHs^CHNHsCuCb is studied 
by neutron scattering in magnetic fields exceeding the critical field of Bose-Einstein condensation 
of magnons. Commensurate long-range order and the associated Goldstone mode are detected and 
found to be similar to those in reference spin-dimer materials. However, for the upper two massive 
magnon branches the observed behavior is totally different, culminating in a drastic collapse of 
excitation bandwidth beyond the transition point. 



Bosc-Einstcin condensation (BEC), such as the super- 
fluid transition in liquid 4 He 1], is the emergence of a 
collective quantum ground state in a system of interact- 
ing Bosons. The condensate is characterized by a macro- 
scopic order parameter that spontaneously breaks a con- 
tinuous U(l) symmetry. For BEC to occur at T > 0, 
the Bosons should be able to freely propagate in 3 di- 
mensions (3D) Q. In one dimension BEC is forbidden 
even at zero temperature. In a striking example of di- 
mensional crossover, even weak 3D coupling can enable 
BEC in ID systems, where the normal state itself results 
from the unique ID topology. A realization of this pe- 
culiar quasi- ID case was proposed only recently and 
involves the condensation of magnetic quasiparticles in 
weakly coupled antiferromagnetic (AF) spin ladders. 

Magnetic BEC can occur in a variety of spin systems 
In gapped quantum magnets, for example, an exter- 
nal magnetic field drives the energy of low- lying magnons 
to zero by virtue of Zeeman effect, prompting them to 
condense at some critical field H c [5(. This transition 
is fully equivalent to conventional BEC. The rotational 
0(2) = U(l) symmetry is spontaneously broken by the 
emerging AF long range order. At T = the density 
of magnons is zero for H < H c , and vanishingly small 
just above the transition. The phenomenon is therefore 
described in the limit of negligible quasiparticle interac- 
tions. To date, such transitions were mainly studied in 
materials composed of coupled structural spin clusters 
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[Tof . [llj ]. The condensing quasiparticles are 
then local triplet excitation that propagate due to inter- 
cluster interactions Il4[ . This is in contrast to the orig- 
inal model of Ref. [3| that deals with coupled extended, 
translationally invariant objects. Their disordered "spin 
liquid" normal state is a direct consequence of ID topol- 
ogy [IH, Even for weak coupling the magnons are 
fully mobile in ID, rather than localized. Is the physics 
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FIG. 1: Time-of- flight neutron spectra measured in IPA- 
CUCI3 in the ID Haldane-gap spin liquid phase (a) and the 
3D ordered BEC phase (b), at H = and H = 11.5 T, re- 
spectively. Solid white arrows indicate the linearly dispersive 
gapless Goldstone mode. 



of the field-induced BEC in quantum AF spin ladders 
any different from that in local-cluster spin systems? 

In the present work we address this issue experimen- 
tally, through a neutron scattering investigation of a 
prototypical spin ladder material (CH 3 )2CHNH 3 CuCl3 
(IPA-CUCI3). This compound almost exactly realizes the 
original theoretical model of Ref. Q. Its spin ladders 
are built of magnetic 5* = 1/2 Cu 2+ ions and run par- 
allel to the a axis of the triclinic PI crystal structure. 
Conveniently, each one can be viewed as a "composite" 
Haldane spin chain [l7j |: pairs of S = 1/2 spins on each 
rung are strongly ferromagnetically (FM) correlated and 
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FIG. 2: (Color online) Background-subtracted neutron spec- 
tra measured in IPA-CuCl3 at the ID AF zone-center h = 0.5 
in various applied magnetic fields (symbols). The field val- 
ues correspond to the ID Haldane-gap spin liquid phase (a,b) 
and the 3D ordered magnon BEC phase (c). Line shapes are 
entirely due to experimental resolution (solid lines). 



act as effective S = 1 objects 17|,ll8j]. Coupling along the 
legs of the ladders is antiferromagnetic (AF), and trans- 
lates into AF interactions between effective spins in the 
"composite" S = 1 chains. Such chains are gapped [l6| 
and are spin liquids with only short-range correlation. 
For IPA-CuCl 3 the energy gap is A = 1.2 meV [13, El]. 

Excitations from this quantum-disordered ground 
state are revealed in inelastic neutron scattering (INS) 
experiments that directly probe the pair spin correlation 
function 5(q, u>). Fig. [Iji shows a time-of-flight (TOF) 
spectrum collected on a 3 g deuterated IPA-CUCI3 single 
crystal sample at T = 100 mK in zero magnetic field us- 
ing the Disc Chopper Spectrometer (DCS) at NCNR and 
6.68 meV fixed incident energy neutrons. The magnon, 
with a steep parabolic dispersion along the a axis and 
gap at the ID AF zone-center h = 0.5, is clearly visi- 
ble. A is small compared to the chain-axis magnon band- 
width, but is considerably larger than transverse band- 
widths along the c (0.4 meV) and b axes (< 0.1 meV) 
[13] ■ Since our main purpose will be to understand the 
special role that the AF spin ladder structure plays in 
the BEC phase of IPA-CUCI3, we shall be comparing our 
results to those found in literature for TICUCI3 [8j, a 
prototypical AF spin-dimer compound. There A is small 
compared to magnon dispersion bandwidths in all 3 di- 
rections [2^, [21]]. To date, TICUCI3 is the only material 
for which the spectrum of spin excitations in the BEC 
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FIG. 3: (Color online) (a) Measured gap energies in IPA- 
CuCb as a function of applied magnetic field, (b) Measured 
fiefd dependencies of the (0.5, 0, 0) magnetic Bragg reflection 
at two different temperatures (symbols). The y-offset cor- 
responds to the actual background for T — 50 mK and is 
arbitrary for the T = 500 mK data. Lines are power-law fits 
over a range of 1 T. Inset: critical exponent f3 as a function of 
the field window used in the least squares fit at T — 50 mK. 



phase has been measured experimentally. 

Figure [2] illustrates the effect of magnetic field on IPA- 
CUCI3. It shows spectra collected at the ID AF zone- 
center h — 0.5 in several fields using the SPINS 3-axis 
spectrometer at NCNR, with 3.7 meV fixed-incident en- 
ergy neutrons, focusing Pyrolitic graphite analyzer and 
a BeO filter after the sample. As the field is turned on, 
the single peak at H = (Fig. [5^) becomes divided into 
three equidistant components (Fig. (2b) . The peak widths 
are resolution-limited. The measured field dependencies 
of the gaps are plotted in Fig. [3]a,, which also includes 
points obtained using the cold neutron 3-axis spectrom- 
eter FLEX at HMI. The gap in the lower mode extrap- 
olates to zero at H c — 9.6 T, where a BEC of magnons 
was previously detected in bulk measurements 22 1. 

As the gap softens, commensurate long-range AF or- 
der sets in and gives rise to new magnetic Bragg reflec- 
tions of type (h + 1/2, k, I), h, k and /-integer. The mag- 
netic structure was determined at H = 12 T in a neutron 
diffraction experiment at the D23 lifting counter diffrac- 
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tometer at ILL, on a 3x2x9 mm 3 single crystal sample, 
using A = 1.276 A neutrons. A good fit to 48 indepen- 
dent magnetic reflections measured at T = 50 mK was 
obtained using a collinear model with spins perpendic- 
ular to the field, aligned parallel to each other on the 
rungs, and antiparallel along the legs of the ladders. The 
refined value of the ordered moment is 0.49(1) /ib- 

The field dependence of the (0.5, —1,0) peak intensity 
measured at T = 50 mK is plotted in Fig. [3b, lower 
curve. To estimate the order parameter critical expo- 
nent (3 we performed power-law fits to the data in a pro- 
gressively shrinking field window (Fig. [3b, insert). The 
extrapolated value is (5 > 0.45, in agreement with ex- 
pectations. Indeed, at T — ► 0, due to vanishing magnon 
density, one should recover the mean field (MF) result 
= 0.5 |,[23|. Any discrepancies between the observed 
and MF behavior become more pronounced at elevated 
T, when magnon density increases, and their interactions 
become relevant. At T = 500 mK, for example, in IPA- 
CuCl 3 we get = 0.25(3) (Fig. [3b, upper curve). BEC 
critical indexes have also been observed under appropri- 



ate conditions in the dimer compound BaCuSi206 11|. 
However, recent work showed that the BEC universality 
of the transition in TICUCI3 is compromised bydevia- 
tions from the Heisenberg model. Anisotropy[l2j and 
magnetoelastic coupling[13j modify the critical indexes 
and account for a small gap in the ordered phase. To 
date, in IPA-CUCI3 we found no evidence of lattice dis- 
tortions at H c or deviations from BEC behavior. 

A key result of this work is a direct measurement of ex- 
citations of the magnetic Bose-Einstein condensate in the 
high-field phase. Data collected above the critical field 
are shown in Fig. [2b (SPINS) and in Fig. [lb (DCS). Three 
distinct excitation branches, two gapped and one gapless, 
are clearly visible, though the overall inelastic intensity is 
reduced compared to lower fields. For each mode the dis- 
persion relations at H = 11.5 T were obtained by fitting 
Gaussian profiles to constant-/! cuts through the data in 
Fig. [lb. The results are plotted in Fig.[4]in solid spheres. 
For comparison, we also plot the dispersion relation of 
the triplet at H — (open circles) [17]. Open squares 
show the dispersion of the middle magnon branch mea- 
sured just below the transition at H = 9 T, on the IN14 
3- axis spectrometer at ILL under similar conditions. 

The gapless mode observed at H > H c at low energies 
(arrows in Fig. [lb) is fully analogous to the phonon in 
conventional BEC, being the Goldstone mode associated 
with the spontaneous breaking of 0(2) symmetry. It is 
a collective excitation of the mag non condensate and has 
a linear dispersion relation [23j. It takes the place of 
the massive quadratically dispersive single-magnon exci- 
tation below H c (Fig. [IJi). At H = 11.5 T the fitted 
velocity of the Goldstone mode vq — 1.74(3) meV is re- 
duced compared to the spin wave velocity [2^] in zero 
field v = 2.9(2) meV 17|. This behavior is qualitatively 
similar to that in TlCuCLj [1, (23[, if one neglects the 
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FIG. 4: (Color online) Dispersion relation of the three excita- 
tion branches measured in IPA-CUCI3 at H — 11.5 T> H c = 
9.7 T (solid symbols). Open circles: magnon dispersion at 
H — 17] Open squares: Dispersion of the middle branch at 
H — 9 T. Lines are guides for the eye. 



tiny anisotropy gap in the latter system [121 ] . which is 
too small to be detected with INS anyway. The sim- 
ilarity also extends to the field dependence of the gap 
energies in those magnon branches that do not soften 
at the transition point. In both compounds the corre- 
sponding slope increases abruptly at H c (Fig. [3b,). A 
bond-operator theoretical treatment of the dimer model 
[23| attributes all these effects to an admixture of the 
higher-energy triplet modes to the condensate. This in- 
terpretation can be qualitatively extended to our case of 
coupled spin ladders. 

While the long-wavelength spectral features in IPA- 
CUCI3 are very similar to those in simple spin-dimer sys- 
tems, the short-wavelength spin dynamics at the zone 
boundary is strikingly different. We find that the two 
massive excitations in IPA-CUCI3 undergo a qualitative 
change upon the BEC transition. As seen in Fig. [lb and 
Fig. H at H = 11.5 T, only 20% above H Cl their band- 
widths are suppressed by over a factor of two. The col- 
lapse occurs abruptly at the critical point: for the middle 
mode there is virtually no change of dispersion between 
H = and H = 9 T (Fig. [4]). Nothing of the sort hap- 
pens in the spin-dimer compound TICUCI3, where the 
bandwidth of the two upper excitation branches evolves 
continuously with field, and is decreased by only 20% at 
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H = 12 T, which is more than twice H c [23[. 

The observed phenomenon can hardly be explained by 
a simple Zeeman shift of quasiparticle energies. Indeed, 
the latter is negligible, as small as ~ 0.2 meV between H c 
and 11.5 T. Instead, we suggest that the abrupt spectrum 
restructuring is is related to the translational invariancc 
of the ground state wave function for an AF spin ladder 
or chain below H c . At H > H c the emergence of long- 
range AF order breaks an additional discrete symmetry 
operation, namely a translation by the structural period 
of the ladder. There is no analogue of this in conventional 
BEC. Depending on inter-cluster interactions that define 
the ordering vector, neither does this necessarily happen 
in spin cluster materials. In particular, in TICUCI3 the 
induced magnetic structure retains the periodicity of the 
underlying crystal lattice. However, in a uniform AF 
spin ladder the extra symmetry breaking is unavoidable, 
regardless of inter-ladder coupling. 

For IPA-CUCI3 the spontaneous doubling of the pe- 
riod implies that at H > H c the wave vectors h = 0, 
h = 0.5 and h = 1 all become equivalent magnetic zone- 
centers. At the same time, h = 0.25 and h = 0.75 emerge 
as the new boundaries of the Brillouin zone. The re- 
sult is a formation of anticrossing gaps for all magnons 
at these wave vectors 24|, where each branch interacts 
with its own replica from an adjacent zone. This trans- 
lates into a reduction of the zone-boundary energy for 
the visible (lower) segments of the two gapped magnons 
in IPA-CUCI3. The additional violated symmetry opera- 
tion is a microscopic one, and therefore plays no role in 
the long-wavelength physics probed at h = 0.5. 

To summarize, any long-wavelength characteristics 
of the field-induced magnetic BEC transition and the 
magnon condensate, such as critical indexes, emergence 
of the Goldstone mode and behavior of gap energies, ap- 
pear to be universal. They are not affected by the ID 
topological nature of the normal state in spin chains and 
ladders, and are very similar to those in local-cluster spin 
systems. In contrast, the short-wavelength properties can 
be significantly different in these two classes of materials. 
In coupled AF spin chains or ladders, unlike in many cou- 
ple dimer systems, and unlike in conventional BEC, the 
transition breaks an additional discrete symmetry. The 
result is a radical modification of the excitation spectrum. 
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